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(54) PROCEDE DE MODE US ATI ON D V UN PROCESSUS PHYSIQUE. 

67) En vue d'automatiser la modelisation d'un processus 
pfiysique a I'aide d'une representation graph ique de type 
schema-bloc, le procede selon i' invention comprend la 
constitution prealable d'une bibliotheque (49) d'objets phy- 
siques symbolises par un bloc et des couples de ports 
d'entree et sortie representant les connexions de Pobjet 
avec son environnement, chaque objet etant associe k plu- 
sieurs modeles mathematiques (50), la determination (41 ) 
des connexions physiques du processus (47) avec son 
Penvironnement, la decomposition (40) hierarchique du 

Crocessus (47) en objets physiques (48) se trouvant en bl- 
liotheque (49), la definition (44) des interconnexions entre 
les ports desdits objets (48), pour chaque objet (48), la se- 
lection d'un modele mathematique (50) en bibliotheque 
(49) suivant la precision recherchee, et la constitution auto- 
matique (45) du modele (51 ) du processus (47) a, partir des 
modeles (50) selectionnes et des connexions entre lesdits 
objets (48). 
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PROCEDE DE MODELISATION D ? UN PROCESSUS PHYSIQUE . 

La presente invention concerne un procede de modelisation 
d'un processus physique a des fins d 1 analyse et de simula- 
tion. 

Elle s 1 applique notamment, mais non exclusivement, a 
l 1 etude du comportement fonctionnel d'un processus physique 
dans son ensemble, a la definition des specifications 
f onctionnelles des organes elementaires intervenant dans le 
processus, de maniere a obtenir un ensemble coherent, ou 
encore, a 1' evaluation des performances du processus, et a 
la definition et la validation des strategies de commande 
du processus . 

A I'heure actuelle, les modeles de processus physiques 
requierent un invest issement important pour leur developpe- 
ment. Par ailleurs, la specificite et la rapide evolution 
des technologies font que ces modeles sont rarement reuti- 
lisables . 




2724744 



On a deja propose des logiciels de simulation et d f analyse 
de processus physiques utilisables dans divers domaines 
d 1 application bases sur l'un ou 1* autre de deux langages 
graphiques principaux, a savoir, le langage schema-bloc 
5 tres utilise par les automaticiens et le langage "bond- 
graph" ou xnultiportes mieux adapts a la description des 
processus physiques , ma is utilisant un forma lisme mo ins 
accessible, 

10 Le langage schema-bloc permet de representer un systdme par 
un assemblage de blocs et de signaux d 1 interconnexions . 
Chaque bloc peut representer n'importe quel type de module 
(localise ou distribue, stabilise ou dynamique) . Cette 
representation est generalement causale, c'est-S-dire que 

15 les differents blocs sont connectes entre eux par des liens 
de cause a effet, et elle d£crit les entrees et les sorties 
de chaque bloc, c»est-a-dire le sens de resolution du 
calcul . 

Ce langage n 1 impose pas de contraintes sur le type des 
20 variables manipulees, un bloc de n entrees et m sorties 
etant caracterise par une relation constitutive explicite 
du type : 

s(t) = f (e(t) ,t) 
e(t) etant les entrees et s(t) les sorties. 

25 

La construction d'un modele graphique par schema-bloc sans 
utiliser une methodologie de developpement s'avere longue 
et coQteuse pour des modeles complexes et presente de 
nombreux inconvenients . 
30 Tout d'abord, les interconnexions entre les differents 
elements n'ont pas de signification physique precise et 
evidente, et ce d f autant plus qu f il existe generalement une 
grande vari^te de variables de couplage, Cet inconvenient a 
pour effet d'empecher d f obtenir des modeles reutilisables. 

35 

Par ailleurs, la structure de tels modeles graphiques etant 
causale, 1 'etablissement d'un modele graphique par schema- 
bloc ne permet pas de mettre en evidence les conflits de 
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causalite, d'oii une difficulty de maltrise du module ainsi 
etabli . 



En outre, lors de la construction d'un modele par schema- 
bloc, le nombre de liens graphiques entre les differents 
elements devient tres vite inextricable, d'otl une faible 
lisibilite et un risque d'erreur important. 

Cependant, cette representation presente I'avantage d'etre 
simple d'utilisation et applicable a tous les domaines. 
Elle est done tres largement utilisee. Elle permet en outre 
de faire apparaltre explicitement toutes les variables 
manipulees quelque soit leur nature* 



Le langage multiportes se base sur la constatation qu*un 
systeme est compose de modules ou de composants dont les 
interactions energetiques ou massiques avec le reste du 
systeme ont lieu a travers un nombre limite de connexions : 
entrees ou sorties de tuyaux, fils electrigues, arbres 
mecaniques, etc . . . 

Selon ce langage, les emplacements des interconnexions sont 
appeles ports. Un systeme ayant plusieurs ports est done 
appele "multiportes" . 

Lorsque deux sous-syst^mes sont relies physiquement par une 
connexion mettant en contact leurs ports, les deux 
variables respectives d 1 effort et de flux des deux sous- 
systdmes sont egales. Ce langage presente I'avantage de 
refleter la structure physique de connexion du systeme et 
d'inclure les lois de conservation au niveau de chaque 
port. 



La presente invention a pour but de cumuler les a vantages 
de ces deux techniques de modSlisation, e'est-a-dire 
d'utiliser les concepts generaux du langage multiportes 
(causality, choix des variables d f etat, conventions de 
signe, ...) en les adaptant, dans un souci de qualite, a la 
representation schema-bloc, qui fait I'objet d'outils 
logiciels de simulation du commerce. 
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A cet effet, ellc propose un procede de mode li sat ion d , un 
processus physique pour simuler et analyser le fonctionne- 
ment dudit processus utilisant une representation graphique 
5 de type schema-bloc, ledit processus faisant intervenir des 
objets, des composants physiques ou des fonctions identi- 
f iables. 

Ce procede est caracterise en ce qu'il comprend : 

10 

- la constitution prealable d'une bibliotheque d f objets 
physiques simples, classes selon les types suivants : 
resistif, capacitif, inertiel, trans forma teur et 
girateur, ou complexes obtenus par association d 1 objets 

15 simples, chaque objet etant symbolise par une representa- 
tion graphique comprenant un bloc et au moins un couple 
de ports comportant un port d 1 entree et un port de sortie 
mis en correspondance, symbolisant les connexions 
physiques de l 1 objet avec son environnement , par lequel 

20 transitent des flux physiques de matiere ou de puissance, 
ou des informations, chaque objet etant en outre associ6 
a au moins un module math^matique de base, simple ou 
complexe caracterise par son degre de precision et sa 
bande de frequences de validite, 

25 

- la determination des connexions physiques par lesquelles 
transitent les flux physiques de matiere ou de puissance, 
ainsi que les informations echanges par le processus avec 
1 • environnement exter ieur , 

30 

- la decomposition hierarchique du processus en sous- 
ensembles physiques, puis en objets physiques se trouvant 
dans la bibliotheque, 

35 - pour chacun des objets physiques composant le processus, 
la determination de ses ports regroupant les variables de 
couplage avec les autres objets physiques const ituant le 
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processus , avec les ports du processus ou du sous- 
ensemble de niveau super ieur, 

- la representation graphique de chaque objet physique 
5 composant le processus par un bloc et des ports, et la 

determination de la causalite de chaque port et de la 
convention de signe de puissance definissant le sens 
positif conventionnel d'ecoulement de la puissance dans 
chaque port de 1' objet, independamment du contexte dans 
10 lequel il s'insere, 

- la definition des interconnexions entre les ports des 
objets composant le processus et la resolution des 
conflits de convention de signe et de causalite, en 

15 tenant compte des conventions de signe et de la causalite 
du processus et celles de certains objets imposees par 
des considerations physiques, 

- la selection dans la bibliotheque de modeles de base 
20 compatibles en causalite, a affecter a chaque objet 

composant le processus en fonction du degre de precision 
a atteindre et de la bande de frequences de validity, de 
mani^re a respecter les connections, les conventions de 
signe et les causa lites precedemment etablies, et 



25 



30 



la constitution automatique du modele mathematique du 
processus a partir des modeles des objets composant le 
processus, et des interconnexions entre les objets prece- 
demment etablies. 



Ce precede permet de realiser un modele graphique hiSrar- 
chise qui constitue une image aussi fidele que possible de 
la composition et de la structure technologique d f un 
processus, dont la modularity facilite devolution et la 
35 maintenance du modele et en ameliore la fiabilit£. 

II permet de reporter la programmation des lois 
const itutives au plus bas niveau, e'est-a-dire au niveau 
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des noeuds terminaux de l'arbre de la decomposition hierar- 
chique, les niveaux superieurs decrivant la structure de la 
decomposition hierarchisee du processus. De cette maniere, 
le concepteur du modele choisit en dernier ressort les 
5 modeles appropries aux objectifs de precision et de bande 
de frequences de validite recherches. II peut ainsi substi- 
tuer un ou plusieurs modeles des objets composant le 
processus par d'autres modules compatibles , sans remettre 
en cause la structure de la decomposition hierarchique et 
10 les connexions entre les objets. 

II permet de reporter au niveau macroscopique le probl&ne 
de la causalite definissant le sens de resolution du 
calcul, c'est-a-dire, parmi toutes les variables identi- 
fiers pour un objet, quel les sont les variables d 1 entree et 
les variables de sortie , ce probleme etant le plus diffi- 
cile a resoudre lors de la modelisation d f un processus* 
En effet, le fait de reporter ce probleme au niveau macro- 
scopique, c'est-a-dire, au niveau du processus lui-meme et 
des sous-ensembles physiques, permet d'une part, de limiter 
le nombre d* objets a prendre en consideration, et d 1 autre 
part, de prendre en compte immediatement la causality 
imposee par certains de ces objets physiques. II s'avere en 
effet que le probleme de la causality d'un processus 
presente une solution pref erentielle ou obliqatoire dans la 
plupart des cas. 

A l^chelle microscopique, le probleme de la causalite est 
deja resolu du fait que l ! on s f en remet a la bibliotheque. 

30 Par ailleurs, ce procede permet de realiser un schema-bloc 
presentant une structure modulaire et f onctionnelle dans 
laquelle les modeles de base representent un ou plusieurs 
objets ou composants physiques. II convient de noter que la 
structure hierarchique globale d'un tel schema est compld- 

35 tement decouplee de la representation locale des objets, 
c f est-a-dire des equations mathematiques decrivant ces 
objets. II en resulte qu f un changement de modele de base 
representant un meme objet physique ne remet pas en cause 



15 



20 



25 
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le modele global representant le processus. Ainsi # la 
precision du modele global peut etre modifiee simplement en 
changeant la precision d'un ou plusieurs modules de base. 

II convient de remarquer egalement que les lois de conser- 
vation de la masse, de la quant ite de mouvement et de 
l^nergie doivent etre verifiees au niveau de chaque sous- 
ensemble, ce qui permet de controler la validite d § un 
modele de base, et done celle du modele du processus. Ce 
procede offre ainsi des facilites de test et de validation 
des modeles ainsi elabores. 



En outre, ce procede permet d* exploiter au mieux le fait 
que ces modeles sont reutilisables en creant une biblio- 
theque comprenant tous les modeles deja realises et que 
I'on complete au fur et a mesure de la construction de 
nouveaux modeles. II permet ainsi au realisateur du modele 
de concentrer son etude sur l 1 analyse globale du processus 
et de s^affranchir de la complexity du modele mathematique 
decrivant le processus, celui-ci etant construit par etapes 
successives a partir de la presentation graphique etablie 
par le realisateur. 



Un mode de realisation du procede selon 1» invention sera 
deer it ci-apres, a titre d f exemple non limitatif, avec 
reference aux dessins annexes dans lesquels : 

La figure 1 represente schema tiquement un calcu- 
lateur mettant en oeuvre le procedS selon 
1 1 invention ; 

La figure 2 est une representation schSmatique de 
l f enchainement des etapes du procede selon 
1 1 invention ; 



La figure 3 represente schematiquement un exemple 
d f un processus que l*on cherche a modeliser ; 
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Les figures 4 et 5 sont des representations 
graphiques du modele du processus obtenues par le 
procede selon 1* invention ; 

Les figures 6 et 7 sont des representations 
5 graphiques des sous-ensembles constituant le 

processus a modeliser ; 

Les figures 8 et 9 montrent deux representations 
differentes d'un meme processus ; 

Les figures 10 et 11 montrent les representations 
10 graphiques, respectivement generique et d6tail- 

lee, d'un sous-ensemble du processus a modeliser. 

Le procede selon 1 1 invention est mis en oeuvre par un 
calculateur du type de celui represents figure 1, compre- 
nant une unite centrale 1 a base de processeurs et de 
15 memoires, des organes de commande tel qu'un clavier 3 et 
une souris 4, un ecran de visualisation 2, et une unite de 
memoire de masse 5 pour stocker notamment une bibliotheque 
comprenant des modeles unitaires et des modeles de proces- 
sus plus complexes deja realises. 

20 

La figure 3 represente un exemple de processus physique que 
l»on peut modeliser gr£ce au procedS selon 1 B invention, Ce 
processus physique est un moteur de type thermique compre- 
nant un papillon des gaz 31, un collecteur d f admission 32, 
25 un moteur thermique 3 3 qui entraine un arbre de transmis- 
sion 35, un dispositif d 1 injection de carburant 34, un 
collecteur d • echappement 36 et un dispositif d % echappement 
37. 

30 Les gaz arrivent en entree du papillon des gaz 31 a une 
pression et une temperature T^. lis ont un debit 

massique mi et un debit enthalpique mi-h^. Les gaz en 
sortie du papillon des gaz 31 sont envoyes au moteur 
thermique 33 par 1 9 intermediaire du collecteur d' admission 
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32 a la sortie duquel agit le dispositif d' injection de 
carburant 34. Lorsque le moteur thermique 33 est alimente 
et en f onctionnement , il applique un couple 38 a I'arbre de 
transmission 35 et emet des gaz d ■ echappement qui traver- 
sent success ivement le collecteur d 1 echappement 36 et le 
dispositif d • echappement 37. Les gaz en sortie du disposi- 
tif d ■ echappement 37 ont une pression P 4 , une temperature 
T 4 un debit massique m 4 et un debit enthalpique m 4 .h 4 . 

Un tel moteur peut etre symbolise grace au proced§ selon 
1 1 invention par un bloc rectangulaire 10 et des f leches 
entrantes lla, 12a, 13a, 14a alignees respect ivement avec des 
f leches sortantes lib, 12b, 13b, 14b s'appliquant sur deux 
bords opposes 15,16 du bloc 10 (figure 4), que l»on peut 
faire apparaitre a I'ecran de visualisation 2 grace a des 
commandes appropriees introduites au moyen des touches du 
clavier 3 ou a l f aide de la souris 4. Chaque couple de 
fleches alignees represente un port du processus a moduli- 
ser . 



Tel que represente figure 2, la premiere etape 41 du 
procede selon 1" invention consiste a definir et a decrire 
ces ports. 

Ces ports peuvent etre de type physique ou inf ormationnel 
et correspondent aux liens du processus avec son environne- 
ment exterieur. Chaque port associant une entree a une 
sortie, est decompose en un couple effort / flux ou pseudo 
e ff°rt / pseudo flux de variables scalaires ou vecto- 
rielles. Les ports physiques peuvent etre de nature 
electrique, mecanique, thermodynamique , hydraulique, etc 
... lis transmettent de l'energie, de la masse ou de la 
quant ite de mouvement, tandis que les ports de type inf or- 
mationnel transmettent de 1 ■ inf ormation et ont une 
puissance negligeable. 



Ainsi, sur la figure 4, le moteur thermique 10 comprend 
deux ports thermodynamiques lla,llb et 12a, 12b de type 
pseudo effort / pseudo flux, representant respect ivement 
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1' admission des gaz de l'exterieur et 1 1 echappement des gaz 
vers l'exterieur. 

Les vecteurs d 1 effort ont pour composantes en entree 11a et 
12a , la pression P^, respect ivement P 4 , la temperature T^, 
5 respect ivement T 4 , et les fractions massiques des consti- 
tuants des gaz X 1X et X 12 , respect ivement X 41 et X 42 - Par 
ailleurs, les vecteurs de flux en sortie lib et 12b ont 
pour composantes le debit massique total m 1# respectivement 
m 4 , le debit massique des constituants des gaz inn et mi2# 
10 respectivement m 4 ^ et m 4 2# et la puissance enthalpique 
m^.h^, respectivement m 4 .h 4 . 

Le moteur thenoique 10 comprend egalement un port mecanique 
13a, 13b associant un couple C mo t en sortie a une vitesse de 
15 rotation w mot en entree, correspondant au transfert de 
puissance vers l*arbre de transmission 5. 

II comprend en outre un port inf onnationnel 14a, 14b compor- 
tant un port d' entree 14a recevant l 1 angle 6 d'ouverture du 
papillon des gaz 31, auquel on a associe un port de sortie 

2 0 fictif 14b pour conserver la symetrie. 

Dans ce schema, l 1 injection de carburant est suppos£e £tre 
une source interne, si bien que le port correspondant est 
interne . 

25 

Au cours de la seconde etape 4 2 du procede, on affecte des 
causalites aux ports 11 a 14- Cette operation consiste a 
choisir parmi les variables associees aux ports, les 
entrees et les sorties du calcul. Or, d'apres le principe 
30 de la causalite en modelisation dynamique, si l'effort est 
en entree, le flux correspondant se trouve en sortie et 
vice versa. 

Cette etape 42 consiste done egalement a associer a chaque 
port, soit une causalite de type effort (flux en entree et 

3 5 effort en sortie) , soit une causalite de type flux (effort 

en entree et flux en sortie) . Cette causalite est symboli- 
s6e sur le schema par des f leches entrantes ou sortantes 
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entre les variables du port et le bloc 10 symbolisant le 
processus. 

Au niveau du processus physique a model iser, la causalite 
des ports est determinee par des considerations physiques* 
Ainsi, le port mecanique 13 est affecte a une causalite de 
type effort, du fait que le moteur thennique est une source 
de couple, modulee par la vitesse w mo t de l'arbre de trans- 
mission 35. II serait en effet un non sens de calculer la 
vitesse a partir du couple. C'est pourquoi la vitesse 
figure en entree du calcul et le couple en sortie. 

En ce qui concerne les ports thermodynamiques 11,12, on 
choisit une causalite de type flux car ils sont relies a 
1 1 atmosphere ambiante qui se comporte comme une source 
d' effort. 



Pour definir completement les ports, il convient ega lenient 
d'affecter arbitrairement une convention de signe a chaque 
port (a !• exception des ports de type inf ormationnel) pour 
indiquer le sens de la puissance, de la masse ou de la 
quant ite de mouvement. 

Une convention positive signifie que les quantites sont 
positives si elles sont revues ou appliquees a l'objet 
modelise, et negatives dans le cas contraire. A l f inverse, 
une convention negative signifie que les quantites sont 
positives si elles sont fournies par l'objet modelise a son 
environnement, et negatives dans le cas contraire. 
Cette convention s' applique aux variables d' effort dans le 
domaine de la mecanique ou la loi de 1' action et de la 
reaction s 1 applique, et aux variables de flux dans les 
autres domaines, oil les lois de continuity s»appliquent 
(convention + si le debit est fourni a I 1 element en 
mecanique des f luides par exemple) . 

Cette convention est etablie sur la base de considerations 
liees a la physique ou a 1' usage. 
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Ainsi, sur la figure 5, on a affecte une convention de 
signe positive (+) au premier port thermodynamique 11 et 
une convention negative (-) au second port thermodynamique 
12 car le flux moyen s'ecoule du premier port vers le 
5 second. 

Une convention de signe negative (-) est appliquee au port 
mecanique 13 du fait que le couple C mot repr^sente le 
couple fourni par le moteur a la transmission. 

10 Une fois que 1 1 environnement du processus a modeliser est 
entiferement defini, il s'agit a la seconde etape 42 du 
procede de resoudre les conflits de causalites. Or, comme a 
ce stade d'avancement de la modelisation, on se trouve en 
presence d'un seul ensemble 10 , aucun probleme de ce type 

15 ne se pose et l f on passe a 1' etape suivante 43. A cette 
etape 43, l f operateur va rechercher si le modele correspon- 
dant a ^ensemble se trouve deja dans une bibliotheque 49 
stockee dans la memoire de masse 5, en fonction du type de 
modele choisi. Si tel n'est pas le cas ou si aucun modele 

20 n'a ete choisi, on passe a l 1 etape 40 qui consiste a decom- 
poser l'ensemble 47 a modeliser, c'est-a-dire, le processus 
physique en sous-ensembles 4 8 en fonction des caracteris- 
tiques physiques du processus. 

25 Ainsi, le moteur deer it precedemment peut se decomposer, en 
se referant a sa structure physique representee figure 3, 
en six sous-ensembles principaux, a savoir : 

- un papillon des gaz 31, 

30 - un collecteur d' admission 32, 

- un bloc moteur 33, 

- un collecteur d 1 echappement 36, 

- un dispositif de detente 37, et 

- un dispositif d 1 injection de carburant 34. 

35 

II convient ensuite de representer chaque sous-ensemble 48 
a l'aide d*un bloc et des ports en procedant de la meme 
maniere que pour le processus physique. C'est pourquoi, le 
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calculateur execute a nouveau les premieres etapes 41 a 43 
pour chaque sous-ensemble 48 que l f on a identifi£. 

Lorsque cela est possible, ou lorsqu'on estime que l^on est 
5 arrive a un niveau de detail suffisant, on peut a cette 
etape selectionner un modele representant le sous-ensemble 
considere. 



Pour cela, on distingue cinq types de modeles elementaires 
10 stockes dans la bibliotheque 49. 

Un premier type R de modele elementaire qui peut etre 
qualifie de resist if, est dissipatif d'energie et est 
associe a aucun etat. II peut etre associe a une causalite 
de type resistance, conductance ou mixte. 

15 

Un second type C de modele elementaire, qualifie de capaci- 
tif, stocke de l f energie potentielle en modifiant ses etats 
qui sont associes a des variables de deplacement. II est 
regi par des equations dif f erentielles dont la resolution 
2 0 correspond a une causality de type integrale, derivee ou 
mixte . 



Un troisieme type I qualifie d'inertiel, emmagasine de 
l'energie cinetique en modifiant ses etats qui sont 
25 associes a des variables de moment. Les modeles de ce type 
sont egalement regis par des equations dif f erentielles dont 
la resolution correspond a une causalite de type integrale, 
derivee ou mixte. 



30 Les modeles elementaires des deux autres types transforment 
l'energie tout en la conservant. Ce sont les transforma- 
teurs TF et les girateurs GY. 

Dans l'exemple precedent, le papillon des gaz 31 peut etre 
35 assimile a un element de type resist if R module par I 1 angle 
d'ouverture e du papillon. 

Le collecteur d' admission 32 est equivalent a un element 
capacitif C. 
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Le bloc moteur 33 est equivalent au niveau macroscopique a 
un transf ormateur TF complexe puisque sa fonction est de 
transformer l f energie thermodynamique en energie mecanique. 
Au niveau microscopique, il peut etre modelise par une 
5 association de modeles de type R, C, I et TF. 

Le collecteur d 1 echappement 3 6 est modelise par un element 
de type capacitif C, tandis que le dispositif de detente 37 
est assimile a un element resistif R. Quant au dispositif 
d 1 injection de carburant 34, il est Equivalent a une source 
10 de debit de carburant liquide ou vapeur. De ce fait, il ne 
comporte qu*un seul port. 

A la premiere etape 41, on definit les ports de chaque 
sous-ensemble 48 precedemment identifie, de la meme manidre 
15 que les ports du processus physique a model iser. Puis, au 
cours de la seconde etape 42, lorsque tous les sous- 
ensembles 48 d'un meme niveau hierarchique ont et£ definis, 
il s*agit d'affecter une causalite a chaque sous-ensemble 
et d'introduire les conventions choisies. 

20 

Cette etape 4 2 revient a real iser une fiche signaletique 
pour chaque sous-ensemble precedemment identifiS en 
association avec 1' etape 43 qui realise 1 1 interconnexion 
des sous-ensembles et qui consiste a interconnecter les 

2 5 symboles representant chaque sous-ensemble 48, en fonction 

de leurs caracteristiques respect ives definies precedem- 
ment. Pour cela, le calculateur 1 cherche a connecter 
chaque port a un autre port associe aux memes variables. 
Cette operation permet d'obtenir le schema-bloc represents 

3 0 figure 6, dans lequel les blocs rectangulaires 61 a 66 

symbolisant les sous-ensembles 48 sont relies les uns aux 
autres par leurs ports de maniere a former une chaine 
presentant certaines analogies avec la representation de la 
f igure 3 . 

35 

Ainsi, sur cette figure 6, le bloc "papillon des gaz" 61 
est connecte au bloc "collecteur d f admission" 62. Le bloc 
"collecteur d 1 admission" 62 est relie a la source 
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"injection" 66 et au bloc "bloc moteur" 63, lequel est 
connecte au bloc "collecteur d 'echappement" 64 ♦ Enfin, le 
bloc "collecteur d' echappement" 64 est connecte au bloc 
"detente" 66. 

La fiche signaletique de chaque sous-ensemble realisee au 
cours de l f etape 42 contient les informations suivantes : 

- le type du sous-ensemble, R, c, 1, TF, GY, ou association 
d» elements simples, 

- la causalite des ports, 

- la convention de signe + ou - affectee au sens d'ecoule- 
ment de la masse, de la quant ite de mouvement ou de la 
puissance, dans chaque port du sous-ensemble, 

- 1' indication si la programmation du modele prevoit ou non 
la reversibilite du flux de matiere dans les ports du 
sous-ensemble , 

- la definition de la base de r§f§rence sur laquelle les 
variables vectorielles doivent etre projetees, ou les 
valeurs de reference dans le cas de variables scalaires. 

La realisation de cette fiche signaletique pour chaque 
ob jet physique permet d' assurer le decouplage entre les 
modeles locaux et le modele global du niveau hiferarchique 
superieur, et d 1 avoir une approche "objet", independante du 
contexte dans lequel le sous-ensemble doit s'inserer. 

L f affectation de la causalite est effectuee par le calcula- 
teur ou l^perateur par une approche hierarchique descen- 
dante qui prend en compte les contraintes imposees par le 
niveau superieur, lesquelles n 1 off rent que peu de degres de 
liberte sur les choix de causalite, en verifiant que les 
causa lites affectees a deux objets a interconnecter sont 
compatibles. 
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Pour cela, le calculateur ou l f operateur commence tout 
d'abord par affecter aux sous-ensembles concernes, les 
causa lites imposees par 1 • environnement du niveau 
super ieur, c'est-a-dire, dans le cas present, du processus 
5 a modeliser. Ainsi, dans ^exemple etudie, la causal ite 
associee aux ports thermodynamiques 11 # 12 de type flux 
implique que la causalite du papillon 31 des gaz et du 
dispositif de detente 37 soient de meme type. Ce qui 
correspond a une causalite de type conductance dans le cas 

10 de modeles elementaires resistifs R. 

Par ailleurs, certains modeles ont une causalite unique qui 
ne depend pas du processus dans lequel ils sont insures. 
C'est le cas d'une source de flux comme le dispositif 
d f injection de carburant 34 qui a une causalite de type 

15 flux. 

L 1 opera teur poursuit si necessaire cette operation en 
affectant ensuite les causalites pref erentielles attachees 
a certains sous-ensembles par prise en compte de considera- 

2 0 tions d' usage liees au metier concerne ou de considerations 

numeriques liees a la representation interne choisie du 
sous-ensemble. Ainsi, dans le cas de modules Elementaires 
capacitifs ou inertiels par exemple, on prefdre en general 
la causalite integrale qui correspond respectivement a une 
25 causalite de type effort ou flux, a la causalite derivee. 

II suffit ensuite de propager les causalites aux autres 
sous-ensembles pour lesquels plusieurs choix sont accep- 
tables, sachant que deux sous-ensembles sont compatibles en 
causalite si l'un est en causalite effort et 1' autre en 

3 0 causalite flux ou inversement. 

Dans le cas o£t 1 1 on aboutit a des conflits de causalite, il 
convient, soit de modifier les causalites prfef Erentielles, 
soit d^jouter des elements de connexion supplementaires 
35 entre les modules de maniere a retablir la coherence. 

Ainsi, tel que represents figure 6, le papillon des gaz 31 
et le dispositif de detente 37 sont en causalite de type 
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conductance ou flux, les collecteurs d 1 admission 32 et 
d 1 echappement 36 en causalite de type integrale ou effort, 
et le bloc moteur 63 en causalite de type flux pour les 
ports thennodynamiques et effort pour le port mecanique. 



tion suivante consiste a affecter a chague port physique de 
chaque sous-ensemble une convention de signe, en tenant 
compte des conventions de signe choisies pour les ports de 
1* ensemble de niveau hierarchique super ieur, c » est-a-dire, 
du processus dans le cas present, et en prenant soin de 
choisir pour deux ports interconnects des conventions de 
signe compatibles. Ainsi, les ports de deux sous-ensembles 
sont compatibles si les conventions affectees respective- 
ment a ces deux ports sont de signes contraires. 
Cette operation permet d'aboutir a la representation de la 
figure 6 ou 7 dans laquelle chaque port des sous-ensembles 
est associe a un signe ( + ) ou (-) . 

Si l f on se trouve en presence de conventions de signe 
incompatibles sur deux ports a connecter, il faut soit 
interposer le signe - sur le fil reliant les efforts dans 
le domaine de la mecanique, ou sur le fil reliant les flux 
dans les autres domaines, soit modifier les conventions de 
signe affectees a ces deux ports* 

II s'agit de determiner ensuite si chaque sous-ensemble est 
a flux reversible ou non, car dans un modele dynamique, il 
peut se produire des inversions de flux locales et on ne 
connait pas tou jours le sens du flux a chaque instant. 

Ainsi, dans l'exemple etudie, seuls les sous-ensembles 
"bloc moteur" 63 et "injection" 66 ne sont pas a flux 
reversible. 

La non reversibilite du sous-ensemble "bloc moteur" 
provient du fait que l'on se trouve en presence d»un objet 
complexe et done que son modele reduit ne s* applique qu'd 
des valeurs moyennes. 



Comme pour la definition des ports du processus, 1' opera- 




Afin d'eviter d'aboutir a des singularites lors de la 
resolution du calcul du modele mathematique global , il est 
en general preferable d'utiliser des modeles a flux rever- 
5 sibles, ce qui implique une certaine redondance dans les 
vecteurs d 1 effort et de flux. 

Lorsque 1 1 on a affaire a des variables vectorielles ou a 
certaines variables scalaires, il est nScessaire de dSfinif 
10 les reperes de reference sur lesquels sont projetes les 
vecteurs , ou de dSfinir des valeurs de reference, comme par 
exemple, dans le cas de tensions electriques ou de 
pressions. 

15 A 1" issue de I'etape 43, on obtient une representation 
graphique dans laquelle les ports des sous-ensembles d'un 
meme niveau hierarchique sont parfaitement definis et sont 
interconnects . 

20 II convient de noter que les liens entre les blocs symbol!- 
sant les sous-ensembles sont tou jours appliques sur les 
bords des blocs paralleles £ la direction genSrale D de 
connexion. En respectant cette regie, on obtient un schSma- 
bloc horizontal comme sur la figure 6, vertical ou croisS 

25 comme sur la figure 7. 

Tel que represents figure 7, la structure graphique de ce 
mode de representation peut refleter la topologie de la 
structure physique du processus a model iser represents 
3 0 figure 3. 



La figure 8 montre un autre exemple d'une representation 
schema-bloc croisee dans laquelle on a assemble le moteur 
thermique 10 tel que decompose figure 6, aux elements 
3 5 mecaniques d'un vehicule a moteur. Ainsi, le "bloc moteur" 
63 assure 1 • entrainement d'une charge 80 constituant le 
vehicule par 1 1 intermediaire d'une transmission 70 modeli- 
see par les elements suivants : 
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l 1 inert ie du moteur 71 , 

un embrayage et une boite de vitesse 72, et 
un arbre de transmission 73. 



Dans ce schema, on a place les sous-ensembles 61 a 66 du 
domaine thermodynamique horizontalement , et les sous- 
ensembles 71,72,73 et 80 du domaine mecanique verticale- 
ment, le sous-ensemble "bloc moteur" 63 assurant l'echange 
entre ces deux parties. Cette representation off re l'avan- 
tage de bien differentier les domaines physiques mis en 
oeuvre par le processus a modeliser. 



En outre, il convient de remarquer sur les figures 6 a 8 
que les croisements de liens sont reduits. Ces croisements 
de liens peuvent etre evites en adoptant une hi6rarchisa- 
tion appropriee. Ainsi, la figure 9 represente les inter- 
connexions entre le moteur thermique 10 et la charge 80 
couplee au moteur gr£ce a la transmission 70. Ce mode de 
representation permet d'aboutir a un diagramme clair et 
synthetique. 



L'etape suivante 44 consiste a ce stade d'avancement de la 
realisation du modele, pour chaque objet ou sous-ensemble 
48 precedemment deer it, soit a selectionner en bibliothdque 
49, un modele mathematique de base en fonction de I'objec- 
tif recherche en termes de precision et de bande de 
frequences, et en fonction de la causal ite precedemment 
etablie, soit a le decomposer en objets physiques plus fins 
(etape 40) . 



A cet effet, la bibliotheque 49 stockee sous la forme d'une 
base de donnees dans l 1 unite de memoire de masse 5 comprend 
les elements suivants pour chaque objet qu f elle contient : 

- le symbole graphique de l 1 objet sous la forme d'un bloc 
et des ports appliques au bloc, definissant les variables 
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correspondant aux ports, la causa lite des ports et la 
convention de signe utilisee, 

- pour chaque modele mathematique associe decrivant 
5 l'objet, la representation graphique sous la forme d'un 
schema-bloc standardise de la structure interne du 
modele, celle-ci etant independante du domaine physique 
concerne , et 

10 - pour chaque modele, la fonction mathematique constitutive 
caracteristique de I'objet represents par le modele sous 
une forme parametrSe dans le cas d'un modele gSnSrique. 

Gr^ce a 1 'application des principes de l'analogie entre les 
15 diffSrents domaines de la physique, la structure interne 
des modeles gSnSriques est independante du milieu physique , 
seules les lois constitutives et la signification des 
variables physiques etant specif iques du domaine concerne. 

20 Cette approche permet d f obtenir un decouplage entre la 
structure generique et le specif ique d'une application qui 
est deer it dans des sous-programmes interchangeables. 

Ainsi, comme represents sur la figure 10, le collecteur 
25 d* admission 32, qui est assimilS a un modele capacitif a 
trois ports 21 a 2 3 de type vector iel en causal ite 
integrale, est represents en bibliotheque de manidre 
generique par un bloc 20 dans lequel les vecteurs de flux 
en entree Fi,F2,F 3 a quatre composantes, sont traites par 
30 quatre sommateurs 25 a 28 qui additionnent les composantes 
de meme rang extraites par des operateurs U(1),U(2), 
U(3),U(4), en tenant compte des conventions de signe affec- 
tees aux ports. Les valeurs delivrees par les sommateurs 25 
a 28 sont ensuite traitees par des integrateurs 1/s, puis, 
35 par la fonction mathematique constitutive $c -1 du modfele, 
avant d'etre envoyees en sortie E par 1 1 intermedia ire un 
multiplexeur MUX. 
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II convient de remarquer que la fonction mathematique 
constitutive peut etre lineaire ou non et est specif ique du 
domaine physique concerne, 

A partir de cette etape 44, il suffit de particulariser le 
modele generique pour obtenir le modele mathematique de 
l'objet physique, comme par exemple celui represents figure 
11* Cette operation consiste a definir la fonction mathema- 
tique constitutive ♦c" 1 modSlisant l f €tat thermodynamique 
(pression, temperature, composition) du gaz dans le volume 
du collecteur d •admission 31, a partir des caract^ristiques 
des flux entrant et sortant, sachant qu'un tel dispositif 
agit a volume constant, c'est-a-dire sans echange de 
travail avec l'exterieur. 



II est a noter que pour deer ire un meme objet, il existe 
une infinite de modeles qui dependent de la precision que 
1 1 on cherche a obtenir et de la bande de frequences que 
l*on cherche a etudier. Parmi toutes les representations 
possibles, on choisira celles qui sont le mieux adaptees a 
la bande de frequences desiree et qui sont compatibles avec 
la causa lite retenue au niveau super ieur. 

Ainsi, par exemple, un bloc moteur peut etre model ise : 

- par une representation du f onctionnement en regime quasi- 
stabilise ; les relations sont alors algebriques et il 
n*y a pas d 1 elements dynamiques, excepte des retards purs 
dQs au transport, ou bien 



- par une representation detaillee microscopique modelisant 
les n cylindres du moteur et le f onctionnement de chaque 
cylindre a l f aide d'un modele complexe compose d'£ laments 
R,C,I,TF,... des domaines thermique, mecanique, thermody- 
namique, . . . 

Si l"on desire af finer davantage la modelisation de l^objet 
en decomposant un ou plusieurs objets 48, on execute a 
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nouveau l 1 etape de decomposition 40 et les etapes 41 a 44 
visant a definir chaque sous-ensemble issu de cette decom- 
position. Bien entendu, cette operation peut etre repetee 
autant de fois que l'on desire. 

Lorsque tous les sous-ensembles 4 8 sont associes £ un 
modele mathematique 50, le procede selon 1 1 invention se 
poursuit par une etape 45 au cours de laquelle Ignite 
centrale 1 elabore le modele mathematique 51 du processus a 
partir des modeles 50 des sous-ensembles 48 et du schema de 
connexion entre les sous-ensembles (figure 6) et des infor- 
mations issues de I 1 etape 4 4 de selection des modeles des 
objets physiques. 

L • etape d f insertion 4 6 des modeles en bibliotheque 49 
consiste a inserer tous les modeles qui ont ete construits 
lors des etapes precedentes, a savoir, les modeles de base 
obtenus en assemblant d'autres modeles de base 50 et le 
modele du processus 51. 

Au prealable, les modeles ainsi realises sont testes et 
valides separement, les lois de conservation de masse 
d'energie et de quantite de mouvement etant applicables au 
niveau de chaque objet. 

Le modele du processus ainsi obtenu peut done etre teste et 
valide facilement, chaque sous-ensemble et chaque niveau 
hierarchique pouvant etre teste et valide separement. Le 
procede selon 1 1 invention permet ainsi de maitriser la 
complexity, ce qui rend plus aisee la realisation de 
modeles complexes fiables. 
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REVEND I CAT IONS 



1. Proc^de de modelisation d'un processus physique 
pour simuler et analyser le fonctionnement dudit processus 
5 en utilisant une representation graph ique de type schema- 
bloc, ledit processus (47) faisant intervenir des objets 
(31 a 37) , des composants physiques ou des fonctions 
ident if iables , 

caracterisg en ce qu'il comprend : 

10 - la constitution prealable d'une bibliothdque (49) 
d' objets physiques simples, classes selon les types 
suivants : resist if, capacitif, inertiel, transf ormateur 
et girateur, ou complexes obtenus par association 
d 'objets simples, chaque objet etant symbolise par une 

15 representation graphique comprenant un bloc (20) et au 
moins un couple de ports (21 a 23) comportant un port 
d»entree (21a a 23a) et un port de sortie (21b a 23b) mis 
en correspondance, symbolisant les connexions physiques 
de 1' objet avec son environnement , par lequel transitent 

20 des flux physiques de matiere ou de puissance, ou des 
informations, chaque objet etant en outre associe Sl au 
moins un modele mathematique de base, simple ou complexe 
caracterise par son degre de precision et sa bande de 
frequences de validite, 

25 - la determination (41) des connexions physiques (11 a 14) 
par lesquelles transitent les flux physiques de matiere 
ou de puissance, ainsi que les informations echanges par 
le processus (47) avec 1 1 environnement exterieur, 

- la decomposition (40) hierarchique du processus (47) en 
30 sous -ensembles physiques (48), puis en objets physiques 

se trouvant dans la bibliotheque (49) , 

- pour chacun des objets (48) physiques composant le 
processus (47), la determination de ses ports regroupant 
les variables de couplage avec les autres objets 

35 physiques (48) constituant le processus (47), avec les 
ports du processus ou du sous-ensemble (48) de niveau 
super ieur, 
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- la representation graphique de chaque objet physique (48) 
composant le processus (47) par un bloc (61 a 66) et des 
ports, et la determination de la causalite de chaque port 
et de la convention de signe de puissance ( (+) ou (-) ) 

5 definissant le sens posit if convent ionnel d • ecoulement de 
la puissance dans chaque port de l f objet (48), indepen^ 
damment du contexte dans lequel il s'insere, 

- la definition (44) des interconnexions entre les ports 
des ob jets (48) composant le processus (47) et la rfesolu- 

10 tion des conflits de convention de signe et de causality, 
en tenant compte des conventions de signe et de la causa- 
lite du processus (47) et celles de certains objets 
imposees par des considerations physiques, 

- la selection dans la bibliotheque (49) de modeles de base 
15 (50) compatibles en causalite, a affecter a chaque objet 

(48) composant le processus (47) , en fonction du degri de 
precision a atteindre et de la bande de frequences de 
validite, de maniere a respecter les connections, les 
conventions de signe et les causal ites prec^demment 

2 0 etablies, et 

- la constitution automat ique (4 5) du modele math^matique 
(51) du processus (47) a partir des modeles de base (50) 
des objets (48) composant le processus, et des intercon- 
nexions entre les objets precedemment etablies. 

25 

2. Procede selon la revendication 1, 
caracterise en ce qu'il comprend 1" entree en bibliotheque 
(49) des modeles de base (50) de chaque objet ou phenomdne 
physique et du modele (51) du processus (47) ainsi determi- 

30 nes. 

3. Procede selon la revendication 1 ou 2 , 
caracterise en ce qu'il comprend la verification et la 
validation du modele mathematique (51) ainsi obtenu, par 

3 5 application des lois de conservation de masse, d'energie et 

de quant ite de mouvement au niveau de chaque objet (48) 
const ituant le processus (47) . 
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4. Procede selon l'une des revendications prece- 

dentes , 

caracterise en ce que la resolution des conflits de causa- 
lite est effectuee par une approche hierarchique descen- 
5 dante, en tenant compte tout d'abord de la causalite de 
l^nsemble de niveau superieur, puis des causalites de 
certains objets (48) imposees par des considerations 
physiques, et enfin des causalites pref erentielles, les 
causalites ainsi dfefinies etant ensuite propagees aux sous* 
10 ensembles pour lesquels plusieurs choix sont acceptables. 

5. Procede selon l'une des revendications pr6c6- 

dentes , 

caracterise en ce que la bibliotheque (49) comprend les 
15 elements suivants pour chaque objet qu'elle contient : 

- le symbole graphique de 1' objet sous la forme d'un bloc 
(20) et des ports (21 a 23) appliques au bloc (20) , 
definissant les variables correspondant aux ports, la 
causalite des ports et la convention de signe utilisee 

20 ((+) ou (-)), 

- au moins un modele mathematique associe a une representa- 
tion graphique sous la forme d'un sch§ma-bloc de la 
structure du modele, et 

- pour chaque modele, la fonction mathematique constitutive 
25 (fcc^ 1 ) modelisant 1 • etat physique de l'objet represents 

par le modele. 



6. Procede selon l'une des revendications prece- 

dentes , 

3 0 caracterise en ce que 1 1 elaboration (43) des modeles mathe- 
matiques (50) des sous-ensembles (48) comprend la selection 
en bibliotheque (49) d'un modele mathematique par objet 
constituant le sous-ensemble en tenant compte de la nature 
des ports du sous-ensemble (48) , de leur causalite, de la 

3 5 convention de signe choisie et du type de modele choisi, et 
la definition de la fonction mathematique constitutive 
( *c _1 ) du modele. 
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7. Procede selon l'une des revendications pr4c€- 

dentes , 

caracterise en ce qu'il comprend, lors de 1 1 interconnexion 
(44) des sous-ensembles (48) , la realisation d'un schema- 
5 bloc dans lequel les symboles de type bloc (61 a 66) repre- 
sentant les sous-ensembles (48) sont interconnects par 
leurs ports de maniere a relier chaque port a un autre port 
associe aux memes variables. 

10 8. Procede selon l'une des revendications prece- 

dentes , 

caracterise en ce que chaque ensemble ou sous-ensemble (48) 
est symbolise par un bloc rectangulaire (10) , chacun de ses 
ports physiques et inf ormationnels etant symbolise par une 
15 fleche entrante (11a a 14a) et une fleche sortante (lib a 
14b) alignees et appliquees respectivement sur deux bords 
opposes (15,16) du bloc (10). 

9. Procede selon la revendication 8, 

20 caracterise en ce que les symboles (61 a 66) de type bloc 
representant les sous-ensembles (48) issus d'une decomposi- 
tion d'un ensemble sont disposes les uns par rapport aux 
autres et relies de maniere a ce que les f leches represen- 
tant les ports soient appliquees sur les bords des blocs 

25 qui sont paralleles a une meme direction (D) . 

10. Procede selon la revendication 8 ou 9, 
caracterise en ce que les symboles (61 a 66) de type bloc 
representant les objets physiques (48) issus d'une decompo- 

3 0 sit ion d'un ensemble sont alignes les uns par rapport aux 
autres 



11. Procede selon la revendication 8 ou 9, 
caracterise en ce que les symboles (61 a 66) de type bloc 
35 representant les objets physiques (48) issus d'une decompo- 
sition d'un ensemble sont disposes les uns par rapport aux 
autres de maniere a refieter la topologie de la structure 
physique dudit ensemble. 
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12. Procede selon l'une des revendications prece- 



dentes, 



caracterise en ce que la definition d'un objet physique 
5 comprend en outre : 

- la definition de son type, simple ( R , C , I , TF , G Y ) ou 
complexe, 

- I 1 indication si son modele mathematique a prevu ou non la 
reversibilite du flux de matiere dans ses ports, 

10 - la definition d'une base de reference sur laquelle les 
variables de couplage de type vectoriel doivent etre 
projetees, ou de valeurs de reference dans le cas de 
variables scalaires . 
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